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Abstract 

Generation of wake fields by a short electromagnetic pulse in a plasma 
with an inhomogeneous background magnetic field and density profile is 
considered, and a wave equation is derived. Transmission and reflection 
coefficients are calculated in a medium with sharp discontinuities. Par- 
ticular attention is focused on examples where the longitudinal part of 
the electromagnetic field is amplified for the transmitted wave. Further- 
more, it is noted that the wake field can propagate out of the plasma and 
thereby provide information about the electron density profile. A method 
for reconstructing the background density profile from a measured wake 
field spectrum is proposed and a numerical example is given. 

1 Introduction 

As is well-known, a short electromagnetic (EM) pulse propagating in an un- 
derdense unmagnetizcd plasma can excite a wake field of plasma oscillations 
(Tajima and Dawson 1979; Gorbunov and Kirsanov 1987). This has interesting 
applications to plasma based particle accelerators (Dawson 1994), photon accel- 
eration (Wilks et al 1989; Mironov et al 1992; Mendonca 2001) and is naturally 
of importance for the general understanding of the interactions between plas- 
mas and radiation. If the plasma is magnetized and the external magnetic field 
non-parallel to the direction of propagation of the exciting pulse, the wake field 
becomes partially electromagnetic and thereby obtains a nonzero group velocity 
(Brodin and Lundberg 1998). 

In the present paper we study wake field generation and propagation in an 
inhomogeneous magnetized plasma. A wave equation for the wake field, includ- 
ing arbitrary inhomogeneities in the particle number density and magnetic field, 
is derived from the cold electron fluid equations and the propagation properties 
are investigated. We address two questions in particular. Firstly, we exam- 
ine the effect of a strong inhomogeneity on the wake field, by introducing a 
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discontinuity in the background magnetic field and density. The longitudinal 
part of the electric field of the transmitted wave can be largely amplified when 
the ratio between the group velocities of the transmitted and incident wave is 
small. The amplification factor for the longitudinal electric field is given and 
analyzed as well as the transmission and reflection coefficients. Secondly, we 
consider to what extent the wake field can propagate out of the plasma, and 
thereby provide information about the background plasma parameters. Since 
the wake field initially has the frequency equal to the local plasma frequency 
also in the magnetized case, this provides a way of extracting information about 
the background electron density profile, i.e. the profile in absence of the wake 
field density oscillations. It turns out that even though wave overtaking - for 
example when a higher frequency part of the wake field passes a lower frequency 
part - may occur, the density profile can still be reconstructed by integrating the 
ray equations of geometric optics backwards. A numerical example is provided, 
where the predicted spectrum of the wake field corresponding to an assumed 
density profile is shown, and a reconstructed profile is calculated. 

The paper is organized as follows: After stating the equations governing the 
wake field in section 2, we derive the wave equation for the longitudinal electric 
field in section 3. The excitation and propagation in a weakly inhomogeneous 
medium are examined in section 4. Then, in section 5, the effects of strong 
inhomogeneities on the wake field, which for example causes field amplification 
and reflection, are studied. Next, in section 6, the spectral properties of an 
electromagnetically generated wake field from a nonuniform density profile are 
investigated. An algorithm for reconstructing the density profile from a mea- 
sured wake field spectrum is given and illustrated with a numerical example. 
Finally, the results are summarized and discussed in section 7. 

2 Basic equations 

We consider a high frequency EM pulse with frequency ujh propagating in a 
cold, inhomogeneous magnetized plasma. We assume the ordering loh ^> oj p , uj c , 
where lu p and uj c = \qHo\/m are the plasma and electron cyclotron frequency 
respectively, q and m are the electron charge and mass, and B = B x is 
the external magnetic field. We let the EM pulse propagate perpendicularly 
to the external magnetic field. The pondcromotive force of the EM pulse will 
generate a "low frequency" wake field mode (which is the low frequency branch 
of the extraordinary mode, or plasma oscillations modified by the magnetic field, 
depending on the choice of terminology) during its path through the plasma. 
The generation mechanism is most efficient if the pulse has a duration of the 
order of the inverse plasma frequency or shorter so that the ion-motion can be 
omitted. In principle, the EM pulse will broaden due to ordinary dispersion, 
decrease its energy and frequency due to the interaction with the wake field, 
etc. These and other effects have been considered in homogeneous plasmas by 
for example Brodin and Lundberg (1998). We will focus on the propagation 
properties of the wake field, however, and for this purpose it turns out that we 
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can forget about the details of the EM pulse. Basically the effect of the EM 
field is to provide a well localized pondcromotive source term in the governing 
equations for the wake field, travelling with almost the speed of light in vacuum. 

The wake field quantities are denoted by index L. We introduce the cor- 
responding vector and scalar potentials Al(z, t) and <f>i,{z,t), using Coulomb 
gauge, and the electron density is written n = no(z, i) + ni,(z, t), where tiq is the 
unperturbed density. Furthermore, the electron fluid velocity v is divided into 
its high- and low frequency part, and we denote the low frequency contributions 
perpendicular and parallel to the direction of propagation with Vii and vl z re- 
spectively. The ponderomotive force of the EM pulse induces longitudinal wake 
field motion, which couple to motion in the y -direction through the Lorentz- 
force, but there is no wake field motion in the direction of the external magnetic 
field, and accordingly we put A £ — A^y and v^i = VL±y. Linearizing in the 
low frequency variables we obtain the following set of equations governing the 
wake field generation and propagation 



Hoqn v L± + [c- 2 d 2 t - dl] A L = (1) 

c~ 2 d z d t (f>L - Hoqno vl z = (2) 

d t v L± - uj c v Lz + — d t A L = (3) 
m 

2 

dt v L z + —d z cj) L + uj c v L ± = gdzlAffl 2 (4) 



in 



d t n L + d z (n v Lz ) = (5) 



3 Derivation of the wave equation 

Before deriving a wave equation for the wake field it is practical to redefine 
*A £ A, ±dA v ± v and S - -(q 2 /2m 2 )d z \A H | 2 . Eqs. §, © and 

(m then reads 



-u 2 p v+[d 2 -c 2 d 2 ]A = (6) 

d t v - Lo c uj- 2 d t ^ + d t A = (7) 

[w- 2 ti$ + l]l> + w c v = S (8) 

Eq. (^) and (0) only give information of how vl z and n l are related to the 

other variables and are therefore omitted at this point. Acting on Eq. ([?]) with 
[df — c 2 d z ] <9 ( -1 and applying Eq. (||) gives 

0=[d 2 -c 2 d 2 +uj 2 p ]v-[d 2 -c 2 d 2 ] M) 

Combining this with Eq. (^|) gives after some rearrangements 

d 2 (ai ! )+^ = j(cj- 1 S) (9) 
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where 



6 = C 2 [10^(1+ LJ p 2 d?)+UJ c U}- 2 ] 

$ - -^ 1 K + 5 t 2 (2 + ^ 2 K 2 +9 t 2 ))] 

7 ^ -[d 2 - c 2 d 2 + . 2 ] 

Eq. (§) provides the starting point for analyzing the propagation properties 
of a wake field generated by an EM pulse under the given circumstances. It 
is somewhat surprising that the evolution of the wake field is governed by a 
single wave equation, in spite that there are two arbitrary background parameter 
functions u p and 10 c . 

In the case of a static and homogeneous background density and magnetic 
field this reduces to 

[df + (u 2 p + u 2 h )d 2 - c 2 d 2 z d 2 - Jl<?d\ + ^] 4, = ^[d 2 - c 2 d 2 + u 2 p ]S (10) 

where io\ = u> 2 + to 2 is the upper hybrid frequency. The left-hand side of Eq. 
(fiXw is the familiar wave operator for the extraordinary electromagnetic mode. 



4 Wake field excitation and propagation 

We will consider wake field propagation in both strongly and weakly inhomoge- 
neous plasmas. By weakly inhomogeneous we mean that the wavelength of the 
wake field is much smaller than the characteristic inhomogeneity length scale. 
Then, to lowest order, derivatives acting on background quantities can be ne- 



glected and the wave equation reduces to Eq. ( |10D with a space and possibly 
a time dependence in tu p and lu c . In a strongly inhomogeneous plasma this 
approximation cannot be applied, and the time evolution of the wake field is 
given by Eq. (||). In this section only weakly inhomogeneous plasmas will be 
considered, and it is illustrative to divide our study of the wake field properties 
into its excitation and its propagation phase. However, we note that the ex- 
citation process considered below will also be of relevance for the next section 
concerning strongly inhomogeneous plasmas, since in that case we will deal with 
wake fields generated in a weakly inhomogeneous plasma that propagates into a 
strongly inhomogeneous region. 



4.1 Excitation 

The excitation of one additional wavelength of the wake field takes place during 
a distance of the order of 2irc/ujp, and - as a basic assumption of ours - the 
variations of uq is negligible on this length scale. Thus as far as the excitation 
process is concerned, the plasma can essentially be treated as homogeneous. The 
solution for the wake field can thus be obtained from previous authors (Brodin 
and Lundberg 1998). Changing to co-moving coordinates £ = z — v g nt, r = t, 
where v g n is the group velocity of the high frequency EM pulse, in Eq. ( |l0| ) and 
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neglecting the small derivatives d 2 and d T d^ and terms proportional to v 2 qH — i 
it reduces to (reinstating the potentials) 



(11) 



This implies 



<t>L = 0Losin[fc p (^ - £ )] 



(12) 



where k p is the wake field wavenumber k p — uj p /v g H, £o is the (constant) position 
of the (short) EM pulse, and c/slo — (qw p /mv g H) A#| 2 d£. The important 
result here, for our purposes, is the determination of the initial value of the wake 
field wave number k p — ui p /v g H, which corresponds to an initial frequency lo p 
(in the laboratory frame). Note that this wake field frequency will in general 
vary with the position of generation. 

4.2 Propagation 

In a weakly inhomogeneous plasma the wavelength of the wake field is much 
smaller than the characteristic inhomogeneity length scale. The wake field prop- 
erties can thus be considered as locally uniform but globally nonuniform, and 
therefore we make the ansatz of geometrical optics (Whitham 1974) 

The local wavenumber and frequency are defined in terms of the eikonal 9{z, t) 
as k = d z 9 and us = —d t 9, respectively and the amplitude i/sq(z, t) is assumed to 
vary slowly with z and t. The local dispersion relation follows then, as a lowest 
order approximation, from Eq. (Mm 



4 2/2 2, 2 ; 2\ , 2 2 ; 2 4 n 

U! - LO (ld h + LO + C k ) + LO h C k + LO = 



(13) 



We note that there is a resonance at or 



LOu 



j h = L0-+L0l 

|[— w c + [loI + 4w 2 ) 1/2 ] and^fl EE \[lo c + (lo 2 c + 4w 2 ) 1/2 ]. The dispersion relation 
has two positive roots. One branch is valid for us > u>r and one for lo^ < lo < lo^. 
As the wake field is generated with the local plasma frequency us p , the wake field 
must belong to the latter branch. Therefore we write the dispersion relation 
from now on as 



W(k,z,t) = { X -yJ; 



— Lo\(?k 2 — LO] 



1/2 



(14) 



where x = ( w l 



J )/2. From the dispersion relation we also derive 



explicit expressions for the wavenumber and group velocity v g 
will use later. They are 



dhus that we 



* = i 

c 



1 1/2 



(15) 
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and 

_ u 2 — ivf c 2 k 

V9 = 2uJ 2 - L0 2 h - UJ 2 - C 2 k 2 ~ ( } 

respectively. 

From the geometric optics approach there follows useful transport equations 
for the wavenumber and frequency. Noting the identity dtk + d z u> = it follows 
that 

%--*W . %=W (IT, 

where dj dt = dt + v g d z . Eqs. ( |l7| ) are referred to as the ray equations. Note 
that for a time-independent medium, the right hand side of the last equation is 
zero, and the wake field propagates with the local wake field group velocity v g 
and with unchanged frequency. 

Since the group velocity may vary along the rays, the energy carried with 
the wake field can be compressed as well as attenuated, and from energy con- 
servation one may expect the field amplitude to vary correspondingly, see e.g. 
Mendonca (2001) for similar effects for ordinary electromagnetic waves in plas- 
mas. The extraordinary mode, however, has several degrees of freedom that the 
energy may vary between, depending on variations in the background parame- 
ters. Therefore the behavior of the amplitudes are not in direct correspondence 
with the variation of the group velocity. For completeness and for future refer- 
ence we state the linear relations between the field variables and ip that follows 
from Eqs. in the weakly inhomogeneous approximation 

4> = -ik- 1 ^ (18) 
u>~ 2 u> c (u) 2 - c 2 k 2 ) / s 



v z — —iojujp 2 ip (20) 

A L = 2 2T2 2^ 21 

Sn = —iknuj~ 2 tp (22) 



5 Reflection and transmission properties 

We now consider the effect of a strongly inhomogeneous region on the wake 
field. We assume that the wake field entering this region was generated in a 
weakly inhomogeneous part of the plasma. Thus the wake field can be taken 
to be uniform when entering the inhomogeneous region, and the variations in 
the wake field frequency u> can be neglected. During these conditions we have 
dfip — ~to 2 4' and, away from the exciting electromagnetic pulse, the wave 
equation @ reduces to an ordinary differential equation 

d 2 z {ail)) +/3ip = (23) 



G 



where 



_ 2-1 —2c 2 2\ 

a = c lu c u p [u h - u ) 

= -a;-V(^ + ^-^ 2 K + ^)) 



Although it is straight forward, at least numerically, to solve Eq. (|23|) for any 
given background density and magnetic field, we simplify the analysis by treating 
the inhomogeneity as a discontinuity, in order to clearly illustrate some of the 
main effects associated with a strong inhomogeneity. We let the discontinuity 
be located at z = 0, and the remaining plasma is assumed to be homogeneous. 
Thus we can make the following ansatz for the wake field 

■02 = 0te i ( fc2 *- wt ) , z > ^ 

where w = uj p \. This ansatz does not apply if we are too close to the exciting 
pulse, or for the fields that were generated in the strongly inhomogeneous re- 
gion, but both these parts of the wake field are assumed to be distant to the 
discontinuity. The subscripts i, r and t stands for the incident, reflected and 
transmitted part, respectively, and the indices 1 and 2 distinguishes quantities 
on the left (z < 0) and right (z > 0) hand side of the discontinuity. 

By integrating Eq. ( p3| ) across z = 0, it follows that aip and d z (aip) are 
continuous over the discontinuity. We define r = ipr/ipi and a = ipt/ipi, and we 
refer to these quantities as the amplification factors for the reflected and trans- 
mitted part of the longitudinal electric field, which essentially are generalized 
Fresnel coefficients. The continuity conditions and the ansatz (E4h imply 



2kl ai (25) 



ki + k 2 ct2 



where ot\ and «2 are the values of a on the left and right hand side of the 
discontinuity, respectively, and 



fci - k 2 
ki + k 2 



(26) 



It follows from Eq. @ that 



S = s (27) 

(lj 2 — c 2 k 2 — uj 2 ) 

is - averaged in time and space - a conserved quantity, i.e. d z S = 0. (Noting 
that S is the time and space averaged z-component of the Poynting vector 
S = i[E x B* + E* x B], this also follows directly from energy conservation. 
Actually from Eq. (^3|) it firstly follows that the conserved quantity is equal to 
ka 2 \ip\ 2 . To see that this is equivalent to Eq. J2~7j ) requires tedious but straight 
forward algebra.) 
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The transmission and reflection coefficients are introduced as T = St/ Si and 
R = S r /Si. Explicitly they read 



T = 



4fcifc 2 



(28) 



(fci + k 2 f 



and 



R = 



(fci - fc 2 ) 2 



(29) 



and they satisfy the energy conservation law T + R = 1 . 

The quantities a, r, T and i? depends on the four parameters ui p i, lo P 2, u> c i 
and lu C 2- We choose oj p \ = lj C 2 = 1 (in normalized units) and, rather than 
presenting complicated surface plots, present one dimensional graphs of the 
dependence of a, r, T and R on uj C 2 for some given values of u> P 2, and vice versa. 
Fig. 1 shows T(o; c2 ), i?(w c2 ) and v g 2 for distinct values of w p2 , and Fig. 2 shows 
T(u> P 2), R(w P 2) and v g 2 for distinct values of w c2 . As can be seen, for most of the 
parameter regime the transmission is close to unity. There are also regions of no 
propagation in Fig. 2, corresponding to parameter values w/ l2 < u) < u)r2- Close 
to these regions the transmission quickly goes to zero and the reflection towards 
unity, because the wave approaches a resonance or a cut-off. The regions to the 
left in Fig. la and 2a corresponds to the branch u) > lur2, i- e - the transmitted 
mode belongs to a different branch of the dispersion relation ( |l3| ) than the 
incident mode. 

In Fig. la, where uj P 2 — 0.5, the left part stretches until w c2 ~ 0.71 where 
the lor cut-off prohibits transmission. At lo C 2 ~ 0.87, T and v g 2 becomes zero 
because of the tou resonance. Fig. lb is the special case for a jump in the 
magnetic field only, i.e. io P 2 = 1. The transmission is everywhere unity although 
the group velocity approaches zero with diminishing w c2 . In Fig. lc, where 
lo P 2 = 2, there is no transmission up to lo C 2 — 3 which is due to the ojl cut-off. 

In Fig. 2a, 2b and 2c the values of co C 2 are 0.5, 1 and 2, respectively. In Fig. 
2a the lor cut-off occurs at lu P 2 ks 0.71, the u>h resonance at lo P 2 ~ 0.87 and the 
u>l cut-off at lo P 2 ~ 1.22. In both Fig 2b and 2c the absence of transmission is 
due to the u>l cut-off. 

The special cases Fig. lb and 2b are particularly interesting, representing a 
jump in the magnetic field only (referred to as case I below) and in the density 
only (referred to as case II below). The transmission remains unity although 
the group velocity goes to zero as lo C 2 — > (Fig. lb) and u P 2 — > (Fig. 2b). 
This means that the energy density entering from region 1 will be dramatically 
amplified in region 2. The question is to which field variables this energy will 
be concentrated. 

Case I. For cj p2 = w pi = 1 the amplification factor for ip is displayed in Fig. 
3a together with v g and v p h- The amplification factor reduces to a — w c1 /cj c2 . It 
should be emphasized that the wavenumber is preserved over the discontinuity 
and thus also the phase velocity, that equals c, is preserved. This property is, 
however, very sensitive to small deviations from exactly constant density. This 
is illustrated in Fig. 3b, where uj P 2 — 0.99, and 3c, where w p2 = 1.01. In 
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these cases the longitudinal field can still be amplified, but not by a large factor 
without also affecting the phase velocity substantially. 

Case II. In the case of uniform magnetic field, lo C 2 — ui c i — 1, the field 
variable that is amplified is the perpendicular electron fluid velocity v±. This 
can be seen by evaluating uxa/vxi using Eq. (Jl^) and the amplification factor 
(p5[). This case does not share the property with case I of preserved phase 
velocity over the discontinuity. In Fig. 2b we have added a small deviation 
to u c so that u C 2 = 1.005 to illustrate that the group velocity can be made 
arbitrarily small. But, in the limit u P 2 — > the group velocity approaches c, as 
required in vacuum. 



6 Density profile reconstruction 

A wake field generated by a short EM pulse in an underdense magnetized plasma 
has the frequency equal to the local plasma frequency, uj p — y / 'q 2 no/eom. Due 
to the presence of the magnetic field it has a nonzero group velocity, and for 
suitable background parameter profiles cut-offs and resonances in the plasma 
are avoided and thus the wake field can propagate out of the plasma. This 
suggests the possibility of gaining information of the density profile no(z) from 
studying the wake field exiting the plasma. 

We assume that the plasma is weakly inhomogeneous so that the results 
in section 4.1 and 4.2 can be applied. Given the ray equations, one may - 
conceptually speaking - treat the wake field as consisting of particles, "photons" 
(dressed photons, that is), created at different times, as the EM pulse plows 
through the plasma. The "photons" have distinct frequencies, depending on the 
position they were "created", and thus also different propagation velocity, v g . 
The equation of motion for the "photons" are Eq. ( |l6| ) and the ray equations, 
Eq. (|l7|). The exact evolution of the wake field, on the other hand, is more 
complicated. This can be understood from that "photons" created at later times 
may be overtaken by "photons" created at earlier times. For our purposes the 
particle-picture is to prefer. 

We assume that the generated wake field spectrum is measured immediately 
outside the plasma boundary. In consistence with the geometric optics approx- 
imation we will treat the measured data as a weakly time dependent spectrum 
with well defined sharp (quasi-monochromatic) peaks. Because of overtaking 
"photons" , the data is not necessarily monochromatic at a given time, multiple 
sharp peaks may occur in the spectrum. Generally, we can express the data as 
a set of distinct frequencies measured at different times, in which case a sharp 
curve can be recognized, see Fig. 4. Due to cut-offs and/or resonances, the 
curve may be discontinuous. 

Given a measured wake field spectrum, as in Fig. 4, the density profile can be 
reconstructed the following way. Discretize the frequency curve into N points 
o/°), a;' 1 ', u;W with corresponding time of detection tw, tw, t^ N '. The 
plasma is discretized into N cells whose positions, z'°), «W and width 

are yet to be determined. The "photon" with frequency u^ N \ detected at time 
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was the last one to exit the plasma. Therefore the plasma frequency in 
the cell at z^ N ' has the value &( N >. Next, we retrace the "photon" lj( n ~ 1 ' ) 
backwards into the plasma to a position z^ N ~^ consistent with the time of 
detection V< N ~ 1 ' and the equation of motions through the already reconstructed 
cell. The plasma frequency in the cell with position z^~ is assigned the value 
^{N— l) rp ne ^ n j rc j "photon" w^- 2 ) is retraced through the cells with plasma 
frequency cjW and w^ -1 ) to a position z^ N ~ 2 ^ and assigns the corresponding 
cell there the plasma frequency uj( n ~ 2 \ This procedure is repeated for all 
"photons" along the frequency curve. 

In order to demonstrate the method we numerically calculate a wake field 
spectrum from an assumed density profile, using the ray equations, see Fig. 5(a), 
and treat this as experimental data from which the density profile can be re- 
constructed. We consider a plasma magnetized such that uj c — 1.1 x w Pi max, 
where tt^max is the maximum value of the plasma frequency. For simplicity we 
normalize such that L = 1, c = 1, where L is the length of the plasma, and let 
the EM pulse enter the plasma at t = Q and exit at t = 1. 

Retracing the spectrum according to the algorithm presented above results 
in a density profile that can be compared with the one we assumed, see Fig. 
5(b). A small numerical error - that can be removed with a finer discretization 
- can be seen. Note that the entire plasma profile cannot be reconstructed. 
The left most points in Fig. 2(b) are missing. This is because the wake field 
generated in this region of low density cannot propagate through the plasma 
since there is a cut-off prohibiting this. The information of this region is already 
missing in the wake field spectrum. 

7 Summary and discussion 

We have considered the propagation of wake fields generated by a short, high 
frequency EM pulse in an inhomogeneous magnetized plasma. A general wave 
equation for the wake field driven by the ponderomotive force of the high fre- 
quency pulse has been derived, Eq. ([)]), and the propagation properties have 
been investigated. If the wake field enters a strongly inhomogeneous region it 
may be largely amplified. The amplification factor for the longitudinal electric 
field and transmission and reflection coefficients have been derived and analyzed 
for a discontinuity in the magnetic field and/or the particle number density. In 
the case of uniform density the amplification factor becomes a = uj c i/uj C 2- This 
special case is also characterized by the fact that there is no reflection and the 
phase velocity of the transmitted wave remains equal to the velocity of light. 
This result may be of relevance for particle and photon accelerators based on 
wake fields. For applications like particle acceleration, it might be desirable 
to finally let the amplified wake field propagate in an unmagnctizcd plasma. 
Simply eliminating the external magnetic field once the wake field is amplified 
affects the frequency and thus the phase velocity in accordance with the ray 
equations ([j"7|). It seems straightforward, however, to match the discontinuity 
and the elimination of the magnetic field so that the phase velocity of the re- 
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suiting field is approximately equal to c, which is the desirable value for particle 
acceleration purposes. 

Furthermore, the spectral properties of a wake field from a plasma with 
nonuniform density have been investigated and a method for reconstructing 
the density profile from a measured wake field spectrum has been proposed 
and illustrated with a numerical example. This result shows that wake fields 
generated by a high frequency EM pulse in principle can be used as a diagnostic 
tool, in magnetized plasmas. The proposed method is based on a mechanism 
substantially different from those of existing techniques, such as interferometry 
and reflectometry (Hartfuss 1998; Hutchinson 1987). It should be possible to 
extend the method to plasmas where the background density varies in time by 
using sequential EM pulses. The most interesting case is that of a strongly 
magnetized plasma (i.e. when the electron cyclotron frequency is larger than 
the plasma frequency), for which almost all of the wake field energy - except that 
generated in a narrow low density region - may propagate out of the plasma. 
The requirement that the length scales of inhomogeneities must be larger than 
the local plasma wave length X p = 2irv g H / <-o P for the results to be valid means 
that it can resolve inhomogeneities of the order 10 8 n~ 1 / 2 m and larger, where n 
is the electron number density in units m~ 3 . 
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9 Figure Captions 



Figure 1. Reflectivity, transmittivity and transmitted group velocity as a func- 
tion of the jump in external magnetic field for oj P 2/ui p \ = 0.5. (Fig. la), for 
u) P 2/uj p i = 1 (Fig. lb) and for uj p2 /uj p i — 2 (Fig. lc), while u> c i/u> p i = 1 in all 
three figures. 

Figure 2. Reflectivity, transmittivity and transmitted group velocity as a 
function of the jump in the unperturbed density for lo C 2/lu c i = 0.5. (Fig. 2a), 
for uj c2 /uj c i = 1 (Fig. 2b) and for w c2 /w c i = 2 (Fig. 2c), while uj c1 /uj p i = 1 in 
all three figures. 

Figure 3. Longitudinal field amplification, phase velocity and group veloc- 
ity of the transmitted wave as a function of the jump in external magnetic field 
for ui P 2/ojpi = 0.99. (Fig. 3a), for lo P 2/oj p \ = 1 (Fig. 3b) and for w p2 /wpi = 1.01 
(Fig. 3c), while uj c1 /lo p i = 1 in all three figures. 

Figure 4. Cartoon picture of a detected wake field spectrum as a function 
of time of detection. The figure also illustrates the disctretization of the curve 
that can be identified from a wake field spectrum. 

Figure 5. Example of a numerically generated wake field spectrum (Fig. 
5a) from an assumed density profile (solid line in Fig. 5b). The reconstructed 
profile is marked with crosses (Fig. 5b). 
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